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Abstract 

This paper is devoted to a long time behavior analysis associated with flow structure interactions 
at subsonic and supersonic velocities. It turns out that an intrinsic component of that analysis is 
the study of attracting sets corresponding to von Karman plate equations with delayed terms and 
without rotational terms. The presence of delay terms in the dynamical system leads to the loss 
of gradient structure while the absence of rotational terms in von Karman plates leads to the loss 
of compactness of the orbits. Both these features make the analysis of long time behavior rather 
subtle rendering the established tools in the theory of PDE dynamical systems not applicable. It 
is our goal to develop methodology that is capable of handling this class of problems. 
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1 Introduction 

The study of von Karman plates in the presence of aerodynamical forces represented by some delayed 
functional is physically motivated [21 These models with delay often arise in the modeling of 

coupled dynamics (e.g., fluid or flow-structure interaction) where the impact of the off-plate dynamics 
can be written as a boundary value of some delayed (flow) potential in the plate equation after a 
sufficiently large time. In fact, this is the case for flow-plate interactions arising in the modeling of 
panels and plates immersed in an inviscid flow (for some discussion, see [B] and |14j). and thus we can 
reduce the study of long-time behavior of solutions to the coupled flow-plate system to the problem of 
a von Karman plate forced by some delay term. 

In the present work, the motivation and significance of studying this class of models derives from 
recent developments in the area of flow-structure interactions, with the goal of attaining good mathe- 
matical understanding of flow-structure dynamics at both subsonic and supersonic flow velocities. It is 
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known from experiment (and also confirmed by numerics), that the potential flow (particularly at the 
supersonic speeds) has the ability of inducing a certain amount of stability in the moving structure. 
This is the case even when the structure itself does not possess mechanical damping mechanisms. If 
one writes down the equations for the interactive system, along with the standard energy balance, this 
dissipative effect is not exhibited at all; quantities are conserved and not dissipated. Thus, there must 
be some "hidden" mechanism which produces this dissipation. It is our goal to shed some light on 
this phenomenon. As it turns out, the decoupling technique introduced in [SJ[S], which reduces the 
analysis of full flow-structure interaction to that of a certain delayed plate model, allows us to observe 
certain stabilizing effects of the flow. These occur in the form of non-conservative forces acting upon 
the structure as the "downwash" of the flow. This idea was already applied to Berger plate models 
[7J [TTJ in the proof of existence of attractors corresponding to the associated reduced plate problem 
with a delayed term. 

In fact, well-poscdness and long-time behavior analyses of nonlinear plate PDEs with delays have 
been treated in [5] (see also [H]): first, in the case of the von Karman model with rotational inertia, 
and secondly, in [Jj [17] , in the case of the Berger model with a small intensity of delayed term (this 
corresponds to a large speed U of the flow of gas - hypersonic). These expositions flesh out the existence 
and properties of global attractors for the general plate with delay in the presence of a 'natural' form 
of interior damping, and then apply this general result to the specific delayed (aeroelastic) force given 
in the full flow-plate coupling. 

It should be noted that the presence of rotational inertia parameter, while drastically improving 
the topological properties of the model, is neither natural nor desirable in the context of flow-structure 
interaction. First, the original model for flow structure interaction describes the interaction between 
the mid-surface of the plate and flow above the plate; and hence should the plate should be treated 
two-dimcnsionally so the equation describing the "downwash" of the flow and the oscillations of the 
plate can agree on the interface. Thus, the rotational inertia term (proportional to the cube of the 
plate thickness, see, e.g., [29]) should be neglected. Secondly, the presence of rotational term changes 
the kinetic energy in the system, forcing a much stronger stabilization mechanism (abstractly, this 
corresponds to the well-known fact that the essential spectrum of an operator can not be moved with a 
compact perturbation). In that case, the stability induced by the flow (a viscous, velocity- proportional 
damping) does not suffice providing a stabilizing effect on the structure. Instead, when the rotational 
inertia term is neglected, the damping secured by the flow alone provides the main mechanism for 
stabilization. This would seem to corroborate physical findings. In view of this, it is paramount to the 
problem at hand to consider the model which does not account for rotational inertia and imposes no 
limiting regimes on the flow velocity parameter U > 0, U ^ 1. 

It is the purpose of this paper to show that the requirement of the rotational term and restrictions 
on the values of U can be eliminated. This is accomplished by resorting to modern analytical tools, 
including compensated compactness methods. These tools, developed over the last few years [14], 
are capable of revealing a mathematical structure of the problem which is consistent with numerical 
and experimental findings. Thus, not only the results, but also the techniques and approach utilized 
in this paper appear to be significant to the discussion of flow-plate interactions and aeroelasticity. 
Accordingly, in this treatment, we focus on the more difficult, and hereto open, non-rotational von 
Karman case, without any restrictions on value of unperturbed speed (0 < U < +oo, U ^ 1). 

The mathematical difficulties which arise in this model force us to consider new long-time behavior 



2 



technologies applied within this framework. In order to provide a glimpse of this, and to demonstrate 
the timeliness of the project undertaken, it suffices to note that the first (and most fundamental) 
difficulty one is faced with is the well-posedness of finite energy solutions corresponding to flow- 
structure interactions - in particularly in the supersonic case. Results and methods found in past 
literature depend critically on the presence of rotational inertial term, see [5j [6] and also |14| Chapters 
6 and 12]. Only recently has this problem has been solved for the non- rotational model for all flow 
velocities [T^] (the subsonic case was also discussed earlier in [HG2I] an d & l so m PS)- Equipped with 
dynamical system structure for the flow-plate interaction, one may then proceed with a study of long 
time behavior. 

The main points of this treatment are: (i) to give a unified approach to the long-time behavior 
of this model, beginning with well-posedness discussions, and working to the existence and finite- 
dimensionality of a global attractor. Our results are novel, and they complete the analysis of long-time 
behavior of von Karman systems with delay by addressing the non-rotational case and the associated 
resulting lack of compactness, (ii) We make use of a recent technique [TTJ [T2] [21 [TS] [21] which allows 
us to obtain the asymptotic compactness property for the dynamical system without making use of 
any gradient type structure of the dynamics (not available in this model, owing to the dispersive flow 
term). This has the added benefit of producing extra regularity of the attractor, and is a less demanding 
approach (the 'traditional' approach of showing quasistability (on the attractor) for von Karman plates 
requires a full characterization of the attractor or assumptions on finitcness of the stationary points 
of the dynamical system and the use of backward-in-time methods, which are not available in systems 
with delay). Lastly, (iii) we make use of new breakthroughs (both in well-posedness and hidden 
compactness) with regard to von Karman flow-plate interactions in order to apply our general results 
on this system with delay to an aeroelasticity model which has received considerable attention from the 
PDE and dynamical systems community, as well as the engineering and aeroscience community, over 
the past 20 years. The results following from the application are also novel, and provide a complete 
treatment of the long-time behavior of a von Karman full flow-plate interaction model. Namely, for 
this model we show the existence of a finite-dimensional compact set in the phase space of the plate 
component. This set attracts all plate trajectories in the case when the corresponding gas flow initial 
data are localized. We emphasize that this result does not assume any damping imposed on the 
system. While stability of flow structure interactions without any damping has been experienced both 
numerically and experimentally, our works appears to be a first rigorous mathematical treatment of 
this phenomena. The key ingredients of the analysis include: exhibit of " hidden" dissipation related 
to the dispersive character of the flow equation along with " hidden" regularity of the boundary traces 
of the flow, the decoupling technique introduced in [6] and applied within the context of recent and 
powerful techniques aimed at showing asymptotic smoothness and quasi-stability |14] without the 
a-priori known compactness or gradient structure. 

We also note that stability and flutter control in flow structure interactions has been also treated 
by several authors at different levels (see, e.g., [2] [3] [21] [23] and the references therein). For instance 
the recent monograph [5] contains a detailed derivation of the flow equation coupled to a linear beam 
equation with the boundary data which involve the Kutta-Jukovsky circulation condition for the flow 
and hence can be applied to elastic wings dynamics. The main trust of the corresponding analysis (see 
also [33J and the literature cited therein) is the characterization of unstable aeroelastic modes arising 
in Possio equation which is linear integral equation describing aerodynamic pressure on the structure. 
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The associated spectral analysis is focused on finding unstable modes in the linear dynamics. The 
obtained results [21 |33] provide specific information on the flutter speed. In contrast, our results 
show that long time behavior of a nonlinear flow-plate model can be reduced to a finite dimensional 
attracting set. More precise information on the stability/unstability of finite dimensional orbits would 
require in depth study of the resulting nonlinear finite dimensional dynamical system (which can be 
chaotic in the supersonic case). 

The paper is organized as follows. In Section[5]we describe the general flow-plate interaction models 
and its relation to the model with delay. We also state well-poscdness results here (Proposition 12. fOl) 
and sketch a proof. This result allows us to define a corresponding evolution semigroup. In Section [3] 
we state and discuss our main result on long-time dynamics of the delayed model (Theorem 13. ip . 
We also show what consequences this result yields for long-time dynamics of the general flow-plate 
interaction model (Theorem 13. 3|) . The next section, Section 2] is central and devoted to the proofs of 
the main results. In Section [5] we briefly discuss the proof of the reduction theorem (Theorem 12.31) 
which rigorously ties the full flow-plate dynamics to the evolution of the von Karman plate with delay. 
Finally, in the Appendix we establish some needed properties of the delayed (aerodynamic type) force. 

Notation: For the remainder of the text we write x for (x,y,z) <G or (x,y) G O C R?^ y )}, as 
dictated by context. Norms || • || arc taken to be L 2 {D) for the domain dictated by context. Inner 
products in L2(R+) are written (■,■), while inner products in L2(R 2 = dM."\_) are written < •, • >. 
Also, H S (D) will denote the Sobolev space of order s, defined on a domain D, and Hq(D) denotes the 
closure of C^°(D) in the H S (D) norm which we denote by || • ||#s(£>) or || • \\ s .d- We make use of the 
standard notation for the trace of functions defined on R5j_, i.e. for (j> £ iJ 1 (]R^), 7[0] = (j>\ z _ Q is the 
trace of <f> on the plane {x : z = 0}. 

2 Motivation and Description of the Model 
2.1 Flow-Plate Interactions 

The model we begin with describes the interaction between a nonlinear plate with a field or flow 
of gas above it. To describe the behavior of the gas we make use of the theory of potential flows 
(see, e.g., (3] 1201 127] and the references therein) which produces a perturbed wave equation for the 
velocity potential of the flow. The oscillatory behavior of the plate is governed by the second order (in 
time) Kirchhoff plate equation with a general nonlinearity. We consider the von Karman nonlinearity, 
which is used in the modeling of the large oscillations of thin, flexible plates - so-called large deflection 
theory. These equations are well known in nonlinear elasticity and constitute a basic model describing 
nonlinear oscillations of a plate accounting for large displacements, see [25] and also P3J fTS] [25] (and 
references therein). 

The gas flow environment we consider is R+ = {(x, y, z) : z > 0}. The plate is immersed in an 
inviscid flow (over body, z < 0) with velocity U ^ 1 in the negative x-direction. (Here we normalize 
U = 1 to be Mach I, i.e. < U < 1 is subsonic and U > 1 is supersonic.) This situation, for instance, 
corresponds to the dynamics of a panel element of an aircraft flying with the speed U, see, e.g., [2T] [22] . 

The plate is modeled by a bounded domain Q C R{( x ^ = {(x, y, z) : z = 0} with smooth boundary 
dil = r and the scalar function u : SI x R + — > R represents the vertical displacement of the plate in 
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the z-direction at the point (x, y) at the moment t. We focus on the plate with clamped boundary 
conditiontQ. 

Accepting von Karman large deflection hypotheses we arrive at the following system: 
u u + A 2 u + ku t + f(u) =p(x,t) in fix(0,T), 

u = d„u = on 91]x(0,T), (2.1) 
^w(O) = u , u t (0) = ui. 

We take the nonlinearity to be von Karman: 

f(u) = -[u,v(u) + F ], (2.2) 
where F$ is a given forcing function and the von Karman bracket [u, v] is given by 

[u, v] = d 2 ud 2 v + d 2 ud 2 v — 2d xy ud 2 y v \ 

and the Airy stress function v{u) is defined by the relation v(u) = u(it, u), where v(u,w) solves the 
following elliptic problem 

A 2 v(u, w) + [it, w] = in il, d„v(u, w) = v(u 7 w) — on (2.3) 

for given u,w £ Hq(CI). 

For the flow component of the model, we make use of linearized potential theory, and we know 
[3j [21] that the (perturbed) flow potential <j> : — > R must satisfy the perturbed wave equation below 
(note that when [7 = this is the standard wave equation): 

\d t + Ud x ) 2 (f> = A4> inR3 x (0,T), 

m = O ; M0) = <h, ( 2 - 4 ) 
d v <t> = d(x,t) onR^ (Xftf)} x(0,T). 

The strong coupling here takes place in the downwash term of the flow potential (the Neumann 
boundary condition) by taking 

d(x,t) = -[(d t + Ud x )u(x)] ■ l n (x), x e R 2 , 

and by taking in (12.11) the aerodynamical pressure of the form 

p(x,t)=p (x) + (d t + Ud x ) 1 [^} with 7 [0] =0| z=o . (2.5) 

Here above ln(x) denotes the indicator function of il in R 2 . This structure of d(x,t) corresponds to 
the case when the part of boundary z — outside of the plate is a the surface of a rigid body. 



1 The clamped boundary conditions are the most physically relevant boundary conditions for the flow-plate model; 
additionally, clamped boundary conditions allow us to avoid certain technical issues in the consideration and streamline 
our exposition. Other possible and relevant plate boundary conditions in this setup include: free, hinged (or simply 
supported), hinged dissipation, and combinations thereof |14l 1291 [30] ■ 
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This gives the fully coupled model: 



'u tt + A 2 u + ku t + f(u) =p + (d t + Ud x )-f[<j>] 
u(0) = u o ; u t (0) = ui, 
u = d v u = 



on dn x (0,T), 
in R 3 + x (0,T), 



in {1 x (0,T), 



(2.6) 



< 



{d t + ud x ) 2 <f> = 

0(0) = O ; ^(o) = 0i, 

d„</> = -[(d t + Ud x )u(x)] -ln(x) 



{(*,!/)} 



x (0,T). 



In this situation, a complete description of well-posedness would require an in depth discussion 
of strong solutions to the (|2.6[) system, including the semigroup formulation and discussion of the 
generator of the dynamics. In addition, these results are not uniform with respect to the parameter 
value U. We refer the reader interested in these details to [H [T31 HH [35] , see also Q3] and the references 
therein. For this treatment, the key point is the well-posedness of weak solutions to (|2.6[) . These weak 
solutions satisfy the variational formulation as defined in |14| . For the purpose of this work we simply 
cite a recently obtained [TB] well-posedness result which attests that the dynamical system generated 
by (|2.6[) is associated to a strongly continuous semigroup on the phase space 



The corresponding result proved in |16j is stated below. 

Theorem 2.1. Let < U ^ 1, k > and F 6 H 3+S (tt), p Q e L 2 (Cl). With reference to the 
system defined in 12. 6}) and any initial data yo = (uo, u%; <f>o, 4>i) G H, there exists a unique solution 
y(-) £ C([0,oo);-ff) which is represented by a strongly continuous semigroup T t : H — > H , y(t) = T t yo, 
t > 0, with the estimate 



where C (R) and uir are positive constants. 

Remark 2.2. When < U < 1 the semigroup T t in Theorem 12.11 is stable in some extended space 
H, i.e. there is a space H 2> H such that ||T f ?/||^ < C(R), t > 0. The above estimate (valid also for 
k = 0) owes its validity to the nonlinear effects [XHl H3] ■ It is not valid for the corresponding linear 
semigroup. 

Various past considerations (se, e.g., [M]) of systems like (|2.6p have made use of an explicit solver 
for the flow. In such an approach, we may rewrite the system above as a von Karman system with delay 
of the form in the earlier sections. Reducing the flow-plate problem to a delayed von Karman plate 
is the primary motivation for this treatment and allows long-time behavior analysis of the flow-plate 
system, which is considerably more difficult otherwise. The exact statement of this reduction is given 
in the following assertion: 

Theorem 2.3. Let the hypotheses of Theorem \2.1\ be in force, and («o, Mi; </>o, 4>i) 6 Hq(£1) x L2(Q) x 
iJ 1 (R^_) x L2(M|]_). Assume that there exists an R such that 0o(x) = 0i(x) = for |x| > R. Then the 



H = £T 2 (O) x L 2 (Q) x H^Rl) x L 2 (R\ 



) 



\T t y ~ TtZo\\ H < C(R)e^ t \\y(0) - z(0)\\ H , V\\y \\ H < R,\\*o\\n < R , 
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there exists a time t^(R, U, CI) > such that for all t > V? the weak solution u(t) to t2.6\) satisfies the 
following equation: 



u tt + A 2 



u + ku t - [u, v(u) + Fq] =po- (dt + Ud x )u - q u {t) 



(2.7) 




q u (t) = — ds d6[M$u](x- (U + sind)s,y - scos6,t - s) 



(2.8) 



Here, u is the extensiorQ of u by outside of ft; Me = sin^9 x + cos9d y and 



t* = in£{t : x({7, 6, s) <£ fl for all x £ ft, 9 £ [0, 2tt], and s > t} 



(2.9) 



with x(£7, 6», s) = (x - {U + sin0)s, y - scos9) C R 2 . 

Thus, after some time, the behavior of the flow can be captured by the aerodynamical pressure 
term p(t) in the form of a reduced delayed forcing. This representation has been used in previous 
considerations of long-time behavior of plates and shallow shells [HJ Section 6.6, pp. 312-334] (and 
the references therein). A rigorous proof of this representation can be found in p~4j pp. 333-334] for 
the rotational case (when we have additional regularity of the plate velocity ut). For the reader's 
convenience, in Section [S] we provide a sketch of the proof, which extends the arguments given in 
|14j for the rotational case. This extension is direct, once the following ingredients are accounted for: 
(1) The full system in (|2.6j) generates strongly continuous semigroup (sec Theorem 12. lj) . (2) the von 
Karman bracket is locally Lipschitz on H (see Lemma f2 . 1 1 1) . and (3) the time derivative of the delayed 
term is bounded on H (inequality (|2.13l) in Proposition 12.41) . 

Theorem 12.31 allows us to suppress the dependence of the problem on the flow variable cj>. Here we 
emphasize that the structure of aerodynamical pressure (|2.5p posited in the hypotheses leads to the 
velocity term —ut on the RHS of (|2.7[) . One can be absorb this term into the damping coefficient k on 
the LHS. However, since we have made no assumptions on the value of k, we may strengthen our result 
for the full reduced flow-plate system by henceforth assuming k = and utilizing the natural damping 
appearing in the structure of the reduced flow pressure, i.e., by moving this term to the RHS. 

As we see below, the reduction method above allows us to study long-time behavior of the dynamical 
system corresponding to (|2.6j) (for sufficiently large times) by reducing the problem to a plate equation 
with delay. The flow state variables {<j), 4>t) manifest themselves in our rewritten system via the delayed 
character of the problem; they appear in the initial data for the delayed component of the plate, namely 
u* | , t , Q ^ . Hence the behavior of both dynamical systems agree for all t > t(R, U, Q). By the dynamical 
systems property for the full system (see Theorem 12. ip , we can propagate forward and simply study 
the long-time behavior of the plate with delay on the interval (<r — t* , a + T] for a > t# and T < oo. 

The following proposition motivates the hypotheses imposed below on the delayed force term in 
the von Karman plate model (|2.7[) . 

Proposition 2.4. Letq u (t) be given by 12. 8\) . Then 




(2.10) 



2 This extension of the solution u(t) is possible owing to the clamped boundary conditions. 
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for any ue L 2 {t-t*,t- : Hl{9)). If u £ L l 2 ° c ([-t*, +oo[; H 2 n i?o)(Q)) we afeo /iawe 



t-t* 



u{t)\\1<It, Vi>0, 



(2.11) 



arte? 



t pt 

\\q U (T)\\ 2 dT < C [t*} 2 \\u(T)\\ 2 2 dT, \/t>0. 



(2.12) 




Moreover if u £ C(-t* , +00; H 2 D H^)(n)), we have that q u {t) £ C 1 (R+; 




For the proof, we refer to Section [5] below. 

Remark 2.5. A priori, when ut is in it is clear from f|2 . 10[) that there is a compactness margin 

and we have the estimate 



However, this is not immediately apparent when u t £ £2^) as ||<7"(i)||(j has no such a priori bound 
from above, as in (|2.10|) . Hence, the critical component which allows us a transition from the 7 > 
case (with damping of the form fc(l — jA)ut) to the 7 = case is the hidden compactness of the 
aforementioned term displayed by (|2.13[) . We note that inequality (|2 . 1 3|) represents a loss of one 
derivative (anisotropic - time derivatives are scaled by two spatial derivatives), versus the loss of two 
derivatives in ([230]) . ([2XTj) . and (f2~T2|) . 

2.2 PDE Description of the Plate Model with the Delay 

Below we utilize a positive parameter < t* < +00 as the time of delay, and accept the commonly used 
(see, e.g., [19] or [35]) notation «*(•) for function on s £ [— 1*,0] of the form s h-» u(t + s). We need this 
because of the delayed character of the problem which requires initial data of the prehistory interval 
[— 1*,0], i.e., need to impose an initial condition of the form u| t6 (_ t » j0 ) = r)(x.,t), where r\ is a given 
function on f2 x [—t*, 0]. We can choose this prehistory data 77 in different classes. In our problem it is 
convenient to deal with Hilbert type structures, and therefore we assume in the further considerations 
that t] e L%{— t*, 0; Hq{CI)). Since we do not assume the continuity of 77 in s G [— 1*,0], wc also need 
to add the (standard) initial conditions of the form u(t = 0) = Uo(x) and d t u(t = 0) = u\{x). 
Again, employing von Karman large deflection hypotheses we arrive at the following system: 



Here f(u) is given by (|2.2[) . The forcing term g(it',£) occurring on the RHS of the plate equation will 




< 



u tt + A 2 u + ku t + f{u) + Lu = po + q(u*,t) in fi X (0,T), 
u = d v u = on dn x (0,T), 
it(0) = uo, «t(0) = Ml, 



(2.14) 



u\ te( _ t * fi) = V eL 2 (-t*,Q;H 2 ({l)). 
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encompass the delayed potential of the gas flow and given by the function q : L 2 (-t*,0; H${n)) xl4 
R, which will be specified below. The scalar k > is our damping coefficient, and represents constant 
viscous damping across the full interior of the plate. The operator L encompasses spatial lower order 
terms which do not have gradient structure (e.g., the term —Uu x in (|2.7[) ). 

Remark 2.6. As it was already mentioned above, the basic plate model we consider may include a 
rotational inertia term (see, e.g. |29j or |f 4j). corresponding to the parameter 7 > and accompanying 
damping parameters k\, &2 > 0. This leads to a plate equation of the form 

(1 - jA)u u + A 2 u + (fci - k 2 A)u t + f(u) + Lu = p Q + q{u\ t). 

Recall that the parameter 7 corresponds to rotational inertia in the filaments of the plate, as discussed 
in the Introduction. These kind of delay models were studied in and [5] and [HI Sections 3.3.1 and 
9.3.1]. We also note that in the case when 

/(«) = /o (/j Vu ( x) ! 2dx ) 

in (|2.14[) , with an appropriate C 1 function f , we arrive to the Berger plate model with delay which 
was studied in [71 [17]. 

Now we formulate our standing hypotheses; we begin with those responsible for well-posedness of 
the model in ([2~14| : 

Assumption 2.7. • We suppose f(u) = — [u, v(u) + Fq], where the functions Fq and po possess 
the properties: 

F (x) G H 3+S (n) for some 6 > 0, p (x) G L a (fi). 

• The linear operator L : i7g(fi) — > £2^) is continuous. 

• t)4 q(v, t) is a continuous linear mapping from Li{—t*, 0; H 2 (£l)) x R+ to L2(£l) possessing the 
property: 

\\q{u\t)\\ 2 < C ( \\u{T)\\ldT, Vt>0, Vu G L l 2 C {[—t* , +00 [; H 2 (fl)). (2.15) 
Jt-f 

Additional hypotheses are needed for long-time dynamics: 
Assumption 2.8. • The linear operator L : Hq~ (Q) — > £2^) * s continuous for some 5 > 0. 

• q(v,t) possesses the (additional) property: 

H^M*,^!! 2 ^ < C [ \\u{T)\\l_ a dT for some < a < 2, (2.16) 

with anyt>0 and u G L l 2 oc {[-t*, +00 [; H 2 (n)). 

• We assume that the generalized time derivative 9t[g(it*,i)] belongs to H~ 2 (Vl) a.s. for any u G 
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C{— t* , T: Hq (fi)) with the following estimate holding for any ip G H^tt): 

\<dt[q{u t ,t)}^>\<c(\\u{t)\\ 2 + \\u(t-t*)\\ 2 + J |Kt + r)|| 2 ^ l^lla-, (2.17) 
for some n > 0. 

Remark 2.9. By Proposition^!]^?/, t) = q u (t) given by (f^T5|) satisfies both Assumptions HTH and |2~51 
Roughly speaking, the conditions in (|2.15|) - (|2.17|) mean that the delay time of the system is distributed 
in the interval [— i*,0], with density which is absolutely continuous with respect to Lebesguc measure. 
This observation also implies that a delay term of the form 

q(v,t) = J Q(t,T)v(r)dT, 

where Q(t, t) is a family of linear operators from Hq (fl) into L 2 (fl) satisfying appropriate hypotheses, 
could be studies from the point of view of this treatment. We also note that we will use the estimate 
in (|2.17p to derive a result on 'hidden' compactness of the term 




< q(u T , t) , u t (r) > dr, 



which is arrived at via integration by parts in time. 



2.3 Well-Posedness of the Plate Model and Energy Relation 

Long-time behavior analysis of the delayed system depends on the well-posedness of suitably de- 
fined weak solutions which generate a dynamical system on the phase space H = H^n) x L 2 (fl) x 
L 2 (-t*,0;H$(n)). 

Well-posedness of weak solutions has been addressed [5] and [HI Section 3.3.1, pp. 189-192; 221- 
222] via the Galerkin method, see also [3 [17] in the case of Berger plates. In what follows we summarize 
and complement relevant results. 

We take a weak solution to (|2.14p on [0, T] to be a function 

u e L^o, T) H^n)) n w£(o, t- l 2 (Q)) n L 2 (-t*, o ; ff 2 (f>)) 

such that the variational relation corresponding to (|2.14[) holds (see, e.g., [131 (4.1.39), p. 211]). We 
now assert: 

Proposition 2.10. Let Assumptions \2.7\ be in force. Then, with initial data 
(u , m.jj) e H = H%(n) x L 2 (Q) x L 2 (-t*, 0; H 2 (ft)), 
problem has a unique weak solution on [0, T] for any T > 0. This solution belongs to the class 

C(0,T;i7 2 (fi))nC 1 (0,T; J L 2 (fi)) 
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and satisfies the energy identity 
S(t) + k [ \\u t (T)\\ 2 dT = g(s)+ I <q(u T ,T),u t (T) > dr+ I < p Q - Lu(t) , w t (r) > dr , (2.18) 



where the full ( not necessarily positive ) energy has the form 

<?(u,u t ) = \{\\u t \\ 2 - < [u,F ],u> } +U,(u) (2.19) 

with 

U4u)^^{\\Au\\ 2 + ^\\Av(u)\\ 2 }. (2.20) 

As stated in P31 Section 4.1.6, p. 221], the proof of Proposition 12.101 requires only minor modifica- 
tions with respect to the proof of the related result in Theorem 3.1.1 [Mj on p. 190. Since Theorem 
3.1.1 deals with rotational models (7 > 0), in order to handle the effect of nonlinear term we rely on 
the sharp regularity of Airy's stress function v(u), given below: 

Lemma 2.11. The function v(u,w) defined in i2.3\) satisfies 

1. |v(«,l0)| W 2,oo(f2) < C||u|| 2 |M|2- 

2. The map u,w — > v(u, w) is locally Lipschitz from Hq(Q) x Hq(Q) —> W 2,oc (fl). 

This lemma easily implies that f(u) in (|2.14[) given by f(u) = — [u, v(u) + Fq] is locally Lipschitz. 

Since the topology for the velocity u t is now L,2{fl), (rather than Hq(FI), as in the rotational case 
7 > 0), we use inequality (|2.15[) rather than (3.3.5) p. 189 [M]. This modification allows us to repeat 
the arguments in P3J , in order to conclude with the statements of Proposition 12.101 For this we also 
make use of the first part of the following lemma: 

Lemma 2.12. We denote q u (t) = g(w',t). Let Assumption \2.7\ be in force. Then 



t r t 

U* I IL./„M|2,J, , , / IL, /,\||2 



<q u (T),u t (T) >dr < Ce~H* / ||u(r)||^r + e / ||u t (r)||^r, Ve > 0, Vt G [0, T], (2.21) 



for any u e L 2 {-t* ,T; H 2 {n)) nWi{Q,T; L 2 (n)). 

If in addition, we assume A s sumption 1 2. 81 then there exists 77* > such that for every e > we 
have the estimate: 

f <q u (T),u t (T)>dT <e f \\u(r)\\ 2 dT + C(t*,e)-(1 + T) sup \\u(t)\\ 2 _^, Vf G [0, T], (2.22) 
Jo J-t* [o,t] 

for any u G L 2 (-t*,T; H 2 (n)) n C*(0, T; H 2 -'i*(n)) n C* 1 ^, T; L 2 (0)). 
For the proof we refer to Section HI 

In order to consider the delayed system as a dynamical system with the phase space H we recall the 
notation: u f (s) = u(t+s), s G [—t*, 0]. With the above notation we introduce the operator St ■ H h-> H 
by the formula 

S t (u , ui,rj) = («(t), ut(t), u*), (2.23) 
where u(t) solves (|2.14p . Proposition 12.101 implies the following conclusion 
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Corollary 2.13. St : H — > H is a strongly continuous semigroup on H. 

Proof. Strong continuity is stated in Proposition 12.101 The semigroup property follows from unique- 
ness. To prove continuity with respect to initial data we use the following assertion. 

Lemma 2.14. Suppose u l (t) for i = 1,2 are weak solutions to {2.1$ with different initial data and 
z — u 1 — u 2 . Additionally assume that 

\\ul(t)\\ 2 + \\Au*(t)\\ 2 <R 2 , i= 1,2 (2.24) 

for some R > and all t £ [0, T]. Then there exists C > and an = an(t*) > such that 

\\z t (t)\\ 2 + l|A2(t)|| 2 < Ce **{||A(uS - u 2 )\\ 2 + \\u\ - u 2 \\ 2 + j ^ \\r,\r) - V 2 (r)\\ 2 2 dr} (2.25) 

for all t e [0,T]. 

Proof. We have that z solves the following problem 

'ztt + A 2 z + kz t + fiu 1 ) - f(u 2 ) = q z (t) - Lz, 

z = d u z = on dfl , (2.26) 
^z(0) = z a E H 2 (n), z t {0) = z x € L 2 (0), z| ( _ t . i0) 6l 2 (-f,0;^(H)), 

where as above we denote q z (t) = q(z t ,t). Let 

E z (t) = ±{\\Az(t)\\ 2 + \\*t(t)\\*}, ( 2 - 2 7) 
then making use of the the energy equality for the difference z 

E z (t) + k f \\z t \\ 2 dT = E Z (Q) ~ I <f{u 1 )-f{u 2 ),z t >dT+ f <q Z {T),Z t (T)>dT 

Jo Jo 

< Lz(t), z t (r) > air 

Jo 

we have the following: 

E z (t) +k f \\u t \\ 2 dT < E z (0) + f [\\fiu 1 ) - f{u 2 )\\ + \\Lz\\] \\z t \\dr\ + 



< q z , z t > dr 

From here, we apply the estimates in (|2.21j) for q z with e = 1/t*, and note the locally Lipschitz 
character of the von Karman nonlinearity f(u), yielding 

E z (t) + kJ^ \\u t \\ 2 d T <C (e z (Q) + J \\z{r)\\ldn) +C(R,t*) J E z (r)dT. 

Hence, Gronwall's inequality yields the desired estimate in (|2 . 25[) . □ 
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Using Lemma r2.14l we obtain that 



\\Styi-S t y 2 \\l i <Ce a * t \\y 1 -y2\\ 2 n 

for Styi = (u l (t), u\(t), [u 1 ]*) with m = (uo,u\,r]), such that (|2.24|) holds. This allows us to conclude 
the proof of Corollary 12. 131 □ 



3 Statement of Main Results 

Our main results deal with (1) long-time dynamics of the system (5t,H) generated by (|2.14[) and (2) 
its connection with the flow-structure dynamics governed by 



Theorem 3.1. Let both Assumptions [S77| and \2.8\ be in force and (St,H) be the dynamical system 
generated weak solutions to the system in (|2.14|) with k > on the space H = Hq(CI) x Lzift) x 
Li{~ t* i 0; Hq (17)) with evolution operator given by (|2.23l) . Then the system (S^H) has a compact 
global attractor A of finite fractal dimension. The attractor can be characterized as the set of all 
bounded full trajectories. Moreover, the set A has additional regularity; namely, any full trajectory 
y(t) = (u(i),w t (i),u*) C A, t e R, has the property that u G L^R; £f 4 (f2) n H$(Q.)) and u t G 
L«,(R; H 2 {n)). 

We recall (see, e.g., [TJ [H [28l [34]) that a global attractor A is a closed, bounded set in H which is 
invariant (i.e., St A = A for all t > 0) and uniformly attracts every bounded set B, i.e. 

lim dn{S t B\A} = 0, where d H {S t B\A} = sup dist H (y, A), (3.1) 

for any bounded BeH, 

The fractal (box-counting) dimension dim/ A of A is defined by 

,. . .. lnn(A,e) 
dim/ A = limsup ■ 



^ ln(l/e) ' 

where n(M, e) is the minimal number of closed balls in H of the radius e which cover the set M. 

Remark 3.2. We note that this type of additional regularity of solutions from the attractor mentioned 
in Theorcm l3.1l is not possible in the case 7 > 0, owing to the fact that the principal term in the equation 
is A 2 it — jAutt- In this case, the presence of this term disallows the use of elliptic regularity theory 
(applied to the biharmonic term) on elements of the attractor. More importantly, when 7 > 0, in 
order to obtain the strong attractiveness property a much stronger damping is necessary In order 
to stabilize the kinetic part of the energy one will have to introduce —jAut as a damping term (see 
section 9.3 in |14j). The point we want to stress is that in our case only u t as a damping term is needed. 
And, in fact, it is this term that will be generated by the flow from "thin air" . As a consequence, the 
plate (in the full flow-plate system) will require no mechanical damping at all. 

Having established the quasicompact character of the delayed potential q u as in (|2.8|) (showing that 
it satisfies the conditions ()2.15p — (|2.17p 'l we can now apply Theorem 13. ll to the von Karman flow-plate 
model in (|2.6p . 
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Theorem 3.3. Suppose < U ^ I, k > 0, F € i7 3+5 (fi) and p e £2^)- T/ierc i/iere exists a 
compact set °^ C Hq(£1) X £2^) = H of finite fractal dimension such that 

Urn dn((u(t),u t (t)),^r) = lim inf (||u(i) - i/ ||^ + ||u t (t) - ^i|| 2 ) = 

t-voo t->oo (i/o.i/ijG* 

/or any weafc solution (u,ut; <f>, (j>t) to \2. 6\) with initial data 

(u^u^fo,^) e H = H^Q) x L 2 (Q) x H 1 ^) x L 2 (R 3 + ) 

which are localized in (i.e., 0o( x ) = </"i( x ) = /or |x| > R for some R > 0). Additionally, we have 
the additional regularity C (H 4 (fl) n H$(Q)) x H 2 (il). 

Proof. The proof of this result follows from rewriting the dynamical system (T t , H) generated by (|2.6p 
as the delayed system in (|2.7p . The latter is possible for sufficiently large times by Theorem 12.31 
Since the delayed potential q u was shown to satisfy Assumptions 12.71 and 12.81 we may apply our main 
result, Theorem 13.11 to the dynamical system generated by the weak solution to (|2.7|) on the space 
H = Hq(SI) x L 2 (Vt) x L 2 (— t*, 0; Hq(Q)). This yields a compact global attractor A C H of finite 
dimension and additional regularity; we then take to be the projection of A on H , which concludes 
the proof as in [14]. □ 

Remark 3.4. We here reiterate that the above result holds in the absence of imposed damping, i.e., 
with k = in (|2.6p . Utilizing the natural damping appearing in the reduced flow pressure, we see that 
in the case of 7 = 0, the flow naturally provides a stabilizing effect to the dynamics in that it yields 
the existence of the compact attractor. This is not the case when 7 > 0, as the nature of the damping 
must be (necessarily) stronger (see Remark l3.2p . 

Remark 3.5. It should also be noted here that because we have rewritten our problem (|2.6p as a 
reduced delayed plate, and additionally changed the state space upon which we are operating, the 
results obtained on long-time behavior will not be invariant with respect to the flow component 
of the model, i.e. our global attractors will be with respect to the state space H, as defined above. 
Again, the data in the form of the delayed term u|/_ t * i0 ) contains the information from the flow itself. 
Obtaining global attractors for the full state space corresponding to (u,ut;4>,4>t) S Hq(£1) x L 2 (^l) x 
i7 1 (M^) x £2(1^) is not a realistic task from the mathematical point of view. There is no damping 
imposed on the system, thus the flow component evolves according to the full half space, unconstrained 
dynamics. The obtained result on the structure (without damping) is the best possible result with 
respect to both the underlying physics and mathematics of the problem. 

Previous Literature and New Challenges: Nonlinear PDEs with delays have been considered 
in various sources (see [3 5) and references therein). In relation to plate equations with delayed aero- 
dynamical type pressure, we note that [14] provides a rather complete analysis of the delayed von 
Karman plate in the presence of rotational terms (and application to flow-plate interactions) , and the 
aforementioned references pj [T7] deal with the plate with delay in the presence of the Berger nonlin- 
earity. In this latter references, it is assumed that the parameter in front of delayed term is suitably 
small. The analysis in the more recent reference |14j also, by and large, applies to the Berger plate (in 
the rotational case). 
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One should stress at the outset that the problem is challenging, even in the rotational case. This is 
due to the fact that the underlying system is intrinsically non-gradient (both the delay term and the 
term Ud x u provide non-conservative and non-dissipative terms that contribute to the loss of gradient 
structure). In view of this, the existence of attractors requires a priori information on a uniformly 
absorbing set. The presence of delay terms along with non-conservative terms makes the latter task 
challenging [7]. 

The references pertaining to plates with delayed terms primarily with the rotational case, i.e. the 
plate equation discussed in Remark l2.6l In particular, for 7, ki,k 2 > 0, well-posedness of weak solutions 
is established via the Galerkin method [14j Sections 3.3.1, 4.1.6]. These solution are shown to generate 
a dynamical system the state space H = Hq(XI) x Hq(Q.) x L,2(— t*, 0; Hq(Q.)). Then, exploiting the 
compactness of the term < q w 1 u t >, i.e. making use of the duality pairing (-ffg, H ), dissipativity 
of the dynamical system can be shown (in much the same way which we utilize below), followed by 
asymptotic smoothness. In this case, however, asymptotic smoothness is arrived at in a straightforward 
way, which additionally exploits the compactness of the von Karman nonlinearity (with respect to the 
energy identity) in the case where u t in Hg(fi). 

The mathematical hurdles arising in the 7 = case (where u t <E L 2 (fi)) begin at the outset 
with well-posedness of weak solutions; indeed, many well-posedness and long-time behavior analysis 
[31)1 [3171 IH)] are dramatically complicated when 7 = 0. Thus, it is no wonder that, to date, the non- 
rotational von Karman plate with delayed terms has not been considered. The reasons for this are 
clear: (i) the methods of studying long-time behavior via an approach making use of the combination 
of Theorems 14.11 and 14 . 2 1 f stated below) is relatively recent. In particular, dealing with the von Karman 
nonlinearity is especially difficult outside of the use of Theorem 14. 2 1 (ii) In addition, the initial studies 
[8] of delayed von Karman plates took place before results on the sharp regularity of the Airy stress 
function were available. These results are critical in this treatment, (iii) Lastly, it is clear by inspection 
that additional properties of the delayed force in (|2.14[) must be accounted for when the rotational 
inertia term is absent. However, owing to a gap in well-posedness results for flow-plate interactions, 
it was unclear what these properties - translated into abstract assumptions - should be. A recent 
observation about the 'hidden compactness' of the reduced delayed potential (derived from an inviscid 
potential flow) yielded insight into what assumptions are reasonable in line with previous analysis of 
von Karman plates with delay. 

In proving finite-dimensionality and smoothness of the attractors, the criticality of the nonlinearity 
and the lack of gradient structure prevents one from using a powerful technique of backward smoothness 
of trajectories [TJ I26[ I14j . where smoothness is propagated forward from the equilibria. Since the 
attractor may have complicated structure, the structure of the attractor is not characterized by the 
equilibria points. In order to cope with this issue, we take the advantage of novel method that is based 
on density and exploits only the compactness of the attractor. 

4 Proof of Main Result 

We first outline the steps utilized to obtain the main result stated above. 

• We begin with recalling the key results on dissipative long-time dynamics for non-gradient sys- 
tems. 
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• We then recall (in the form adapted to the delay case) cite the primary estimates which have 
been used in previous long-time behavior considerations for von Karman plates in the past. 
Noting that we cannot make use of the 7 > approach to long-time behavior (owing to the loss 
of compactness of the term < q u (t),ut(t) >). We begin by exploiting a different assumption 
('hidden' compactness of this term when integrated in t in the energy relation ([2TB]) ). 

• We then use a similar, modified Lyapunov functional as that in [21 p. 480] on the dynamical 
system to show that it is dissipative. 

• After obtaining the necessary compactness estimates, we synthesize them to produce a pointwise 
energy estimate which allows us to make use the abstract Theorem 14.21 to obtain asymptotic 
smoothness of the dynamical system associated with weak solutions to (|2.14|) . At this point, we 
make use of abstract Theorem 14. II to conclude that the dynamical system possesses a compact 
attractor in the space H. 

• In the last step, we revisit our estimation in the asymptotic smoothness section to obtain the so- 
called quasistability estimate on the attractor utilizing its compactness; this allows us to apply 
Theorem 14.31 The application of this theorem gives the finite dimensionality and additional 
smoothness of the attractor. 

4.1 Preliminaries on Dissipative Dynamical Systems 

We recall notions and results from the theory of dynamical systems (see, e.g., [T] 19] 128] 134]). 

One says that a dynamical system (S^H) is asymptotically smooth if for any bounded, forward 
invariant set D there exists a compact set K C D such that 

lim d H {S t D\K} = 

t— y+00 

holds. An asymptotically smooth dynamical system should be thought of as one which possesses local 
attractors, i.e. for a given forward invariant set Br of diameter R in the space H there exists a compact 
attracting set in the closure of Br, however, this set need not be uniform with respect to R. 

A closed set B C H is said to be absorbing for (S t ,H) if for any bounded set D C H there exists 
a to(D) such that StD C B for all t > to. If the dynamical system (St,H) has a bounded absorbing 
set it is said to be dissipative. 

In the context of this paper we will use a few keys theorems (which we now formally state) to prove 
the existence of a finite dimensional global attractor. First, we address the existence of attractors and 
characterize the attracting set: 

Theorem 4.1. Any dissipative and asymptotically smooth dynamical system (St, H) in a Banach space 
H possesses a unique compact global attractor A. This attractor is a connected set and can be described 
as a set of all bounded full trajectories. 

For the proof, see [1] or [34] . 

Secondly, we state a useful criterion (inspired by [26] and proven in [12], see also [HI Chapter 
7]) which reduces asymptotic smoothness to finding a suitable functional on the state space with a 
compensated compactness condition: 
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Theorem 4.2. Let (St,H) be a dynamical system. Assume that for any bounded positively invariant 
set BcH and for all e > there exists a T = T € ^b such that 

\\Stxi - StX2\\h < e + ^ t ,B.T{x\,x 2 ), x t E B 

with $ a functional defined on B x B depending on e, T, and B such that 

liminf liminf ^e,T,B{x m ,x n ) = 

m n 

for every sequence {x n } C B. Then (St,H) is an asymptotically smooth dynamical system. 

In order to establish both smoothness of the attractor and finite dimensionality, a stronger estimate 
on the difference of two flows is needed. We now cite [14] pp. 381-387] and also [15]. Note that we 
have used a specialization of the cited theorem which utilizes the special structure of the state space 
in the problem at hand: H = H 2 (£l) x x L 2 (— t*, 0; Hq(Q,))); the theorems cited above are more 

general: 

Theorem 4.3. Let A be a global attractor for (St, IV). 

x\,x 2 G B C H where B is a forward invariant set for the flows StXi. Assume that the following 
inequality holds for all t > with positive constants Cx, C 2 , to 

\\S t Xi - S*z 2 ||h < Cxe-^Wxx - x 2 \\ 2 H + C 2 max || Ul ( T ) - u 2 (r)\\ 2 Hi (4.1) 

re[0,t] 

for any x\,x 2 G A, where StXi = (ui(t) , dtUi(t) , u\) and Hq(Q) C i?i C L 2 (Q) is compactly embedded. 
Then the attractor A posesses the following properties: 

(a) The fractal dimension of A is finite. 

(b) For any x £ A one has dt(Stx) 6 L^R, H). 

The estimate in (|4.1j) is often referred to (in practice) as a "quasistability" estimate. It reflects the 
fact that the flow can be stabilized exponentially to a compact set. Alternatively, we might say that 
the flow is exponentially stable, modulo a compact perturbation (lower order terms). The quadratic 
nature of the lower order terms is important for the validity of Theorem 14.31 

4.2 Technical Preliminaries 

In this section we derive and cite certain energy and multiplier estimates, as well as estimates on the 
von Karman nonlinearity, which will be necessary in the proof of Theorem 13.11 below. 

The following theorem is a case specialization found in [14] Section 1.4, pp. 38-45; Section 9.4, 
pp. 496-497]. These bounds elucidate the local Lipschitz (quasi-Lipschitz) character of the von Karman 
nonlinearity are relatively recent and critical to our nonlinear analysis. 

Theorem 4.4. Let u 1 G @ R (H§(Q)), i = 1,2, and z = u 1 - u 2 . Then for f{u) = ~[u,v(u) + F ] we 
have 

||jV)-/(" 2 )ll-a < C s (l + \\u 1 \\ 2 + \\u 2 \\ 2 )\\z\\ 2 _s <C(S,R)\\z\\ 2 - S for all tfe[0,l]. (4.2) 
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If we further assume that u l £ C(s, t; H 2 (fl)) D C 1 (s, t; L2(£l)), then we have that 

< fin 1 ) f{u% zt >= \j t Q(z) + \P{Z) 



where 



and 



Q{z)^<v{u 1 ) + v{u 2 ),[z,z] > -\\Aviu 1 + u 2 , z)\f 



P{z) = -< ul, [u\v{z)\ >-<u 2 t , {u 2 ,v{z)\ >-<u\ + u 2 tl [z, v^u 1 + u 2 , z)} > . (4.3) 



Moreover, 



<f(u 1 (T))~f(u 2 (T)),z t (r) >dr < C{R) sup \\z\\%_ v + \ [ P{z)dr (4.4) 

T£[s,t] A Js 

for some < i] < 1/2 provided u'(t) <E &r(H$(Q)) for all r E [s,i\. 
The above bounds rely on the equation 

/(u 1 ) - f(u 2 ) = [z, v(u) + F ] + [u 2 , viu 1 ) - v(u 2 )} 

and on the so-called 'sharp' regularity of the Airy stress function v(u) (see Lemma \2. lip . 

We will now make use of the above estimates in producing energy type estimates. 
First, we multiply (|2.14[) by u and integrate over the set f2 x (s, t), making use of clamped boundary 
conditions. This produces the following identity: 



<u t (r),u(r)>\ t + f (\\Au(r)\\ 2 -\\u t (r)\\ 2 ) 



dr 



( < J(w(r)),w(r) > + < q u (r) + p - Lu,u(t) > W 



where as above we use the notation q u (r) = q(u T , r). By standard splitting and interpolation, we 
arrive at 

J (||A W (r)|| 2 - || Ut (r)|| 2 )dT < ej \\u\\ 2 dT + e\\p + 0\\ 2 - J < /(«(t)),u(t) > dr 

+ C [* \\q u (T)\\ 2 _ a dT + C e (l + \t-s\) sup |Kr)||t„ 



reU.t] 



+ I < u t (t),u(t) > I + I < u t (s),u(s) > | 

for all e > and for some i] > 0. This estimate, coupled with the estimates in (|2.16[) and (|4.2[) yield 
the following estimates: 

Lemma 4.5. Let u 1 6 C*(0, T; H 2 (fl)) (1^(0, T; L 2 (Q)) DL 2 (-t*, T; H 2 (fl)) solve $1%) with clamped 
boundary conditions and appropriate initial conditions on [0,T] for i = 1,2. Then the following 
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estimate holds for all e > ; for some r\ > 0, and < t < T : 



f (\\Au\\ 2 -\\u t \\ 2 )dT< e f \\u\\ 2 2 dT + e + C f \\u(r)\\ldT + C(e,t*,T) sup ||u(t)||§_„ 

Jo JO J -t* re\0.t] 



< f(u), u> dr + \ < u t (t), u(t) > | + | < u t (s), u(s) > |. 



Moreover, in the case where we are considering the difference z = u 1 —u 2 of solutions solving \2. 26\) with 
u % (t) G 3§r(H 2 (Q)) for all t G [0, T], we may utilize the estimates in Theorem \4-4\ (which eliminates 
the stand-alone e ) arrive at 

f {\\Az\\ 2 -\\z t \\ 2 )dr< e f\\z\\ldT + C f \\z{T)\\l_ dT + C{e t T,R) sup \\z(r)\\ 2 2 _ v 

Js Js Js-t* re[0,t] 

+ E z {t) + E z (s), (4.5) 

where E z (t) is given by (f2T27j) . i.e., E z {t) = ^{\\Az(t)\\ 2 + ||z f (i)|| 2 }. 

The final class of estimates we need are energy estimates for the z term defines as the solution 
to (|2.26p . Energy estimates for single solutions (making use of the nonlinear potential energy) can 
be derived straightforwardly from (|2.18p . The energy estimate on z, along with the estimate in (|4.5p 
above, will be used in showing asymptotic smoothness for the system. 

The energy relation on [s, t] for z in (|2.26p is given by 

E z (t) + k f \\z t \\ 2 dT = E z (s) - f <f(u 1 )~f(u 2 ),z t >dT+f <q z (T),z t (T)> dr 

J S J S J S 



s: 



< Lz(r),z t (r) > dr 



From this, making use of splitting and Sobolev inequalities, we quickly have for < s < t < T, some 
rj > 0, and all e > 0: 

E z (t) + k f \\z t \\ 2 dT<E z (s)+C(e,T) sup \\z\H_,, + e [ {\\z\\ 2 + \\z t \\ 2 )d T 

Js T£[s,t] Js 

< fiu 1 ) - f(u 2 ),z t > dr + j < <f(T),z t (r) > dr 
In the case when k > by (I4.5|) this implies that 

\E z {t)+CQ I E z dr < E z (s) + C(T,R) sup \\z\\ 2 2 _ n + C f \\z(r)\\ 2 2 ^dT 

^ Js T6[s,i] Js-t* 

- / <f(u 1 )-f(u 2 ),z t >dr+ f <q z t (r),z(r)>dT 

J S J S 

+ \<q z (t),z{t) > | + | <q'(s),z(s) > |. 
Here above we usee the integration by parts formula for the integral with delayed term. Therefore 
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there exist a, > and C(T,R) > such that 

E z (t)+ I E z dT<a (E z (s)+ f \\z{T)\\l_ a d^] +C(T,R) sup \\z\\ 2 2 _^ 

Js \ J s-t* J re[s,t] 



- a x J < f{u x ) - f(u 2 ), z t > dr. 
Taking t = T and integrating over s in [0, T] we arrive at E| 

TE Z (T)+ [ ds [ E z dr < a (e z (0) + / \\z(T)\\ 2 2 _ a dT^) +C{T,R) sup ||«||__ 

Jo V J-t* / TG[0,T] 

- ai / ds / < /(w 1 ) - ,/(u 2 ), z t >dr -a 2 ( < f{u v ) - f(u 2 ), z t > dr. 



Since 



we arrive to the following assertion. 



ds / E z dr > 



T 



E z dr for T > 2t* 



T-t* 



Lemma 4.6. Let u l £ C(0, T; # 2 (ft)) n C^O, T; L 2 (fi)) nL 2 (-t* , T; i^(fi)) sotoe fO^P wi/i damped 
boundary conditions and appropriate initial conditions on [0,T] for i = 1,2, T > 2t*. Additionally, 
assume u l (t) g 3§ji(H 2 (Q)) for all t 6 [0, T]. TTien £/ie following estimates 



E Z {T)- 



T-t* 



E z (r)dT < a o [E z (0)+ \\z(r)\\idT ) + C(T,R) sup ||z||^ (4.6) 

re[0,T] 



ai / ds / < fiu 1 ) - f(u 2 ),z t > dr -a 2 < f{u x ) - f(u 2 ),z t > dr 



hold with ai independent of T and R. 



4.3 Dissipative Dynamical System 

Our next task in order to to make use of Theorem 14. II is to show dissipativity of the dynamical system 
(5*4, H), namely that there exists a bounded, forward invariant, absorbing set. To show this, similar to 
the consideration in |14l Theorem 9.3.4, p. 480], we consider the Lyapunov type function (with S{u, Ut) 
as in (|2~19| and with II* (u) given by (|2~20| ) 

V(S t y) =S(u(t),u t {t))- < q u (t),u(t) > +v{ < u u u > +^\\u\\ 2 ) (4.7) 

+ fx( [ n*(u(s))ds + / ds [ n*(it(r))dr), 

^ J t-t* JO Jt-s ' 

3 With other constants a; and C(T,R). 
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where Sty = y(t) = (u(i), Ut(t), it*) for t > and /n, f are some positive numbers to be specified below. 
In view of the results for the von Karman plate in [14j Section 4.1.1], we have that 

coE(u,u t ) -c<V(S t y) <c x E{u,u t )+ fiCt* J II*(u(t + r))dr + c (4.8) 

for j/ > small enough, where co,ci,c, C > are constant. Here we make use of the notation: 
E(u,u t ) = ^\\u t \\ 2 +n*(u). 

To obtain the above bound, we make direct use of our assumption on the L2 bound in (|2.16[) on 
the term < q u ,u >. Additionally, we here need (and below) a critical lower bound on the potential 
energy, which can be found in various forms throughout [141 p. 49 and p. 132]: 

Proposition 4.7. For any u G H 2 (H) n Hq (f2) we have for any 6 > and any < 77 < 2 

IMlL,<<5(|M| 2 + ||A«( u )|| 2 ) + C 5 , 

where v(u) = v(u,u) is the Airy stress function as defined in (|2 . 3[) . 

In what follows below, we will often make use of the above theorem to give 

< [u,F ],u> +|M| 2 < <5II»(u) + C s , Fo - 

We now compute —V(S t y): 
dt 

^-V{S t y) = ^7-£{t) + v < u tt ,u > +^||u t || 2 + kv <u t ,u> 
dt dt 

- < q u ,u t > - < q£,u > 
d 



|jT n*{u{s)ds + J n*(u(r))drds| 



= — fc||?it|| 2 + < po — Lu, Ut > — < qf, u > 

+ v <u tt + ku t ,u > +^||M t || 2 + |un*(u(t)) - fjH(u(t - t*)) 

+ (it*IL*(u(t)) - n J n*(«(t - r))dr 

Now, we make use of the relation 

u tt + ku t = -A 2 u+p Q + q u (t) + [u, v(u) + F ] - Lu 
owing to (|2.14[) . Substituting this back into the relation above and simplifying yields: 

d r° 

jV(S t y) =-(k- v)\\u t \\ 2 - v\\Au\\ 2 - v\\Av(u)\\ 2 - ^(u(t - <*)) -p] n.(u(i + r))dr 

+ v < [u, Fq],u > +v < q u , u > + < po — Lu, ut > +v < po — Lu, u > 
-<#,u>+Mi* + l)n.(u(i)). 



21 



In the estimate that follows we make use of (a) standard splitting via Young's inequality, (b) the bound 
in Proposition 2171 (c) the assumption that < v < min{l, k}. Then we have for all e > 

-nii(u(t-t*))-n J n*(u(t + T))dT 

+ e|M| 2 + £ (I|A U || 2 + \\Av(u)\\ 2 ) + C e , s , Po , Fo + e\\q»\\ 2 + | < q?,u > \. 

Now, using ([235]) : 



\\q u (t)\\ 2 <cf ||A W (r)||^r, 
Jt-f 



and (|2~T7|) with i 



) = u: 



\<q?(t),u(t) > | < 



\\Au(t)\\ 2 + \\Au{t - t*)\\ 2 + / \\Au(t + T)\\ 2 dT 



+ Ce,t*\Ht)\\ 



we have (again, making use of Proposition ^. 7\i the following inequality for all e > 0: 

jV(S t y) <(-k + v + e)\\u t \\ 2 + (^I±^2 + 2e v)\\Au\\ 2 + (^i±£l + e - „) \\Av(u) 
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- (fj, - 2e)n*(i - t*) -(/j,- 4e) J II, 



(u(t + T))dT + C e ||u(t)|| J + Ce.^.Fo. 



And, for < v < fc, and for /i and e sufficiently small, we have the following lemma: 

Lemma 4.8. For all k > there exist n, v > and c(/i, v,t* 1 k) > and C(/U, i^poi ^b) > smc/i £/iai 

d r° 

J t V ( S tV) < -c{\\u t \\ 2 + \\Au\\ 2 + \\Av(u)\\ 2 +U r (u{t - t*)) + J H,(u(i + r))dr} + C. 

From this lemma and the upper bound in (|4.8[t . we have for some (3 > sufficiently small (again, 
depending on /i and v): 

jV(S t y) + f3V(S t y)<C, t > 0, (4.9) 
The estimate above in (|4.9|) implies (by a version of Gronwall's inequality) that 

F(5 t y)<y(y)e-^ + ^(l- e -^). 

Hence, the set 

y G H : ^(y) < 1 + ^ 
is a bounded forward invariant absorbing set. This gives that (H, St) is dissipative. 
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4.4 Asymptotic Smoothness 

Recall that our dynamical system is (St, H), where St is the evolution operator corresponding to plate 
solutions to (|2.14j) and H = Hq(CI) x L 2 (fi) x L 2 (— t*, 0; Hq (fi)). To show asymptotic smoothness of 
this dynamical system, we will make use of an abstract Theorem 14. 2 1 

To make use of this theorem, we will consider our functional ^ to be comprised of lower order 
terms (compact with respect to H) and quasicompact (/ < f(u r ) — f(u 2 ),z t > dr) terms. On the 
LHS of the above estimate, we see that we need to produce an estimate which bounds trajectories in 
H, i.e. ||(u(t), Ut(t), m*)||h (taking the metric d to be || • ||h)- Such an estimate will be produced below 
by combining our energy estimates produced earlier: 

Lemma 4.9. Suppose z = u 1 — u 2 is as in 112.26)) . with y l (t) = (u z (t) , u t (t) 1 ,u t ' 1 ) andy l (t) £ ^_r(H) 
for all t > 0. Also, let rj > and E z (t) be defined as in \2.21ty . Then for every < e < 1 there exists 
T = T £ (R) such that the following estimate 



E Z (T)+ f \\z{T)\\ldT<e + ^^ T Av\y 2 ) 

JT-t* 



holds, where 



*e,TM(y\y 2 )=C(R,T) SUP \\Z(T)\\1_ V + a, [ <f(u 1 (T))-f(u 2 (T)),Z t (T)>dT 



"2 



re[0,T] 
T 

< f(u\r)) - f(u 2 ( T )),z t (T) > drds 

Proof. It follows from (|4.6|) by dividing by T and taking T large enough. □ 

In Lemma 14.91 above, we have obtained the necessary estimate for asymptotic smoothness; it now 
suffices to show that as defined above, has the compensated compactness condition described in 
Theorem [4~2l 

Before proceeding, let us introduce some notation which will be used throughout the remainder of 
this section and in the following section. We will write 

l.o.t.= sup \\z{r)\\l_ v , F{z) = fiu 1 )- f(u 2 ). (4.10) 

t£[0,T] 

Theorem 4.10. The dynamical system (S t , H) generated by weak solutions to ^2.14\ l is asymptotically 
smooth. 



Proof. In line with the discussion above, we aim to make use of Theorem 14.21 To do so, it suffices to 
show the compensated compactness condition for f^j which we now write as "3/, with e,T, and R 
fixed along with the other constants given by the equation. Let B be a bounded, positively invariant 
set in H, and let {y n } C B C ^^(H). We would like to show that 

liminf liminf ^(y n ,y m ) = 0. 

m n 

More specifically, for any initial data Uq = (uj, w}, 7/ 1 ), U 2 = (u^,u\,r] 2 ) € B (where rf belongs to 
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L 2 (-t*,0 : Hg{Q))) we define 



T 



T r T 



^([/ \[/ 2 ) = / <T(z)(r),z t (T)>dT + / / <F(z(r)),Zt(T) > drds 

JO JO Js 

where the function z = u 1 — u 2 has initial data Uq — Uq and solves (|2.26[) . The key to compensated 
compactness is the following representation for the bracket [HI PP- 598-599]: 

< F(z){r),z t {r) >= ~{ WAviu 1 )]] 2 - \\Av(u 2 )\\ 2 + 2 < [z,z],F > } 
- < [v(u 2 ),u 2 ],ul > - < [^(u 1 ),^ 1 ],^ 2 > . 

Integrating the above expression in time and evaluating on the difference of two solutions z n ' rn = 
w n — w rn with initial data Wq — W™, where vf — 1 w, yields: 

lim lim / T <^(z"' m )(r),z^ m (r)>dT = i{||Ai;(u ) ( S ))|| 2 -||At-(u;(T))|| 2 } (4.11) 

n— >oo rn— >oo J s 2 

lim lim / { < [v(w n ),w n ],wr > + < [v(w m ),w m ],w? > }, 



where we have used (i) the weak convergence in H 2 (£l) of z n ' m to 0, and (ii) compactness of Av(-) 
from H 2 (Q) — > £2^) as in Lemma [2.111 The iterated limit in (|4.11[) is handled via iterated weak 
convergence, as follows: 

lim lim / { < [v{w n ),w n ],wr > + < [v{w m ), w m ], < > } 

n^-oo m— >oo J 

= 2 J* < [v(w),w],vh >= l\\Av(w)(s)\\ 2 - ±\\Av(w)(T)\\ 2 . 
This yields the desired conclusion, that 

r-T 



lim lim I {T < z n - m (T)), 2?' m (r) > dr = 0. 



The second integral term in "J 2 is handled similarly. Since the term l.o.t. above is compact (below 
energy level) via the Sobolev embedding, as a consequence we obtain 

liminfliminf4 f (y n , y m ) = 0. 

m— >oo n— >oo 

This concludes the proof of the asymptotic smoothness via Theorem 14.21 □ 



Having shown the asymptotic smoothness property, we can now conclude by Theorem 14.11 that 
there exists a compact global attractor A C H for the dynamical system (St, H). 
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4.5 Quasistability Estimate 

In this section we refine our methods in the asymptotic smoothness calculation and work on trajectories 
from the attractor, whose existence has been established in the previous sections. 
Analyzing (|4.6p . we may also write 



where J-(z) is given in (|4.10[) . We note that c does not depend on T > min{l, 2i*}, and l.o.t. is of 
quadratic order. 

In order to prove the quasistability estimate (as in (|4.1[0 . we have to handle the non-compact term 
< J 7 {z),z t >. We recall the relation (@3) in Theorem 0~4] if u l G C(s, t; H 2 (n)) n C\s, t; L 2 (fL)) with 
u 1 {t) e ^ R (L 2 (Cl)) for t e [s,t], then 



for some < rj < 1/2. Here P(z) is given by (14.31) . 

Let 7 u i = {(u 1 (t),ul(t), [u 1 ]*) : t 6 1} and 7ll 2 = {{u 2 (t), u 2 (t), [u 2 ] 1 ) : t 6 R} be trajectories from 
the attractor A. It is clear that for the pair u 1 ^) and u 2 (t) satisfy the hypotheses of the estimate 
in (|4.13[) for every interval [s,i\. Our main goal is to handle the second term on the right hand 
side of (|4.13|) which is of critical regularity. To accomplish this we shall use the already established 
compactness of the attractor in the state space H = Hq(Q.) x £2^) x L 2 (—t*, 0; H 2 (Vl)). 

Since for every t£I, the element u\(t) belongs to a compact set in L 2 (£Y), by density of Hq(Q) 
in L 2 (ft) we can assume, without a loss of generality, that for every e > there exists a finite set 
{<pj} C Hq(£1) , j = 1, 2, n(e), such that for all r S R we can find indices ji(r) and j 2 {r) so that 




(4.12) 




(4.13) 



KV) - Jl(T) || + \\u 2 (t) - ^ 2(t) || < e for all ret. 



Let P(z) be given by (|4.3[) with the pair u : (t) and w 2 



(t) and 



(0 3 i, [u 1 , - (fe, [u 2 ,v(z,z)}) - {<j) n +<f) h ,{z,v(u' 1 + u 2 ,z)]) 



where z(t) 



u\t) 



u 2 {t). It can be easily shown that for all ji,j 2 < n(e) 



\P(z(t)) - P jl(T)tMT) (z(r))\\ < eC(A)||z(r)|| 2 



(4.14) 



uniformly in r € R. 
Starting with the 



estimate (1.4.17) page 41 [T4] . 



IIKHII-2 <c|H| 2 _^|Hli+/s, V/3e[o,i) 
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and exploiting elliptic regularity one obtains 

1 1 [it, V(Z, w)}\\- 2 < C||lt||2-/3||[2,li;]||_2 < C||lt||2-/3||2||2-/3 1 ||w||n-/3 1 , ( 4 -15) 

where above inequality holds for any j3, pi £ [0, 1) 

Recalling the additional smoothness of <fij £ Hq({1), along with the estimate in (|4.15|) applied with 
j3 = j3i = r], and accounting the structure of Pj terms one obtains: 

\\P juh (z)\\< C(A)(ll&J|2 + \\4> h \\2)\\z(T)\\ 2 2 _ v 

for some < rj < 1. So we have 

sup||P^,i 2 (^)|| < C{e)\\z{T)\\l_ rj for some < rj < 1, (4.16) 
ji -h 

where C(e) — > oo when e — > 0. Taking into account (|4. 14[) and (|4.16[) in (|4. 13[) wc obtain 

f <F(z),z t > <C(e,T,A) sup ||z(r)||t.„ + e / \\z(r)\\ldr (4.17) 

Js re[s,t] Js 

for all s £ R with rj > and t > s. Considering (|4.17[) and taking T sufficiently large, we have from 

(031 

E Z (T)+ [ \\z(T)\\ldT< 1 (E z (0)+ f \\z(r)\\ldr)+C(A,T,k,t*) sup ||z(r)||t, 

JT-f J-t* r£[0,T] 

with 7 < 1. Thus by the standard argument (see [13]) we finally conclude that for y(t) = (z(t), z t (t), z 1 ) 
[|w(*)IIh < C7(a-, A)||y(0)]|^e- rt + C sup \\z(t)\\1_. 

r£[0,t] 

Hence on the strength of Theorem 14. 3[ applied with B = A and 

H = H%{tt) x L 3 (fi) x L 2 (-t*,0;H$(n)), 

where H\ = H 2 ~ r, (fl), we conclude that A has a finite fractal dimension. 
Additionally, Theorem 14.31 guarantees that 



\\u tt {t)\\ 2 + \\u t (t)\\ 2 < C for all t e M. 
Since u t £ H 2 (tt) C C(f2), elliptic regularity theory for 

A 2 u = —u u — ku t — /(it) — Lu + q(ii*, t) 
with the clamped boundary conditions give that 

\\u{t)\\ 2 4 < C for all t S M. 
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Thus, we can conclude additional regularity of the trajectories from the attractor A C H stated in 
Theorem 13. II 

We have now completed the proof of Theorem 13.11 

4.6 Extensions and Open Problems 

In this section, we briefly mention possible extensions of our results and open problems pertaining to 
the models discussed above. 

• Other homogeneous boundary conditions: hinged, simply supported, free-type, or combinations 



• Nonlinear internal damping, i.e. considering k(u t ) in the plate equation, where £;(•) is a Nemitsky 
operator. 

• Boundary damping via (nonlinear) hinged dissipation (|31j). 

• Other physical nonlinearitics; replacing the von Karman nonlinearity in the considerations above 
with Berger or Kirchoff-type nonlinearity (as discussed in [!|] Q2] ) ■ This should not present 
major technical issues. 

• Convergence to equilibria; one may conjecture that individual trajectories of the full flow-plate 
system converge to single elements of the attractor. However, Dowcll's conjecture (supported 
by numerics) states that only in the subsonic case solutions stabilize to stationary points, and 
in the supersonic case, locally stable periodic (or even chaotic) orbits are possible. Hence, it is 
likely that one can discuss this convergence only in the subsonic case. The principal issue here 
is finiteness of the dissipation integral 



We can easily guarantee this if we have additional plate damping in the coupled system. In this 
case we can achieve stabilization in the same way as it done in |14| for the rotational case (see 
also [TO])- So the issue becomes how to obtain some form of finiteness of dissipation integral in 
the case where full energy of coupled system is preserved (the reduction result Theorem [33] is not 
employed). This issue remains open (see also the corresponding remark in |14[ Section 12.4.2]). 



thereof. 
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5 Flow- Plate Interactions 

In this section we provide sketch of the proof of Theorem 12.31 which is needed for our principal 
application of Theorem 13.11 to the long-time dynamics of the fully coupled model: 

' u t t + A 2 " + ku t + f(u) =p + (d t + Ud x )-y[4>] 
u(0) = u o ; u t (0) = ui, 
u = d v u = 

' (d t + ud x ) 2 4> = Acb 
0(0) = O ; <k(o) = <h, 

d v <j> = -[(8t + Ud x )u(xj\ -ln(x) 

Proof. By Theorem 12.11 the system (|5.1|) generates a strongly continuous semigroup T t on H . This 
is to say that (4>(t),4>t(t),u(t),Ut(t)) = T t (4>o,(j)i,uo,ui),t > 0. The main idea behind the proof of 
Theorem 12.31 (see [21 Section 6.6.5]) is to split gas flow variable <fi in two components: 0(x, t) = 
0*(x, t) + 0**(x,t), where 0*(x, t) solves (|2.4j) with d(x, t) = and 0**(x, t) is solution to non- 
homogenous problem (|2.4j) with the zero initial data 4>o = and 4>\ = 0. Here we note that with 
(0o, 4>i) £ ,ff 1 ( R +) x ^2(R+) one obtains [3! [16] ((/>* (t) , $ (t)) G H 1 (M.\) x L 2 (M+)- Thus, by Theo- 
rem [5j] w e also have that 

Note that £/iis Zas£ property is not valid for a flow solution with Li boundary Neumann datcQ. However, 
the improved regularity is due to the interaction with the plate and the resulting cancelations on the 
interface. Moreover, we also obtain a meaningful "hidden trace regularity" for the aeroelastic potential 
on the boundary of the structure |16) : 

(d t + Ud x H<P] G L 2 (0, T; H-^ 2 (n)) (5.2) 

where T is arbitrary. 

Then, using the Kirchhoff type representation for the solution <fi*(x,t) in (see, e.g., Theo- 
rem 6.6.12 in |14j). we can conclude that if the initial data 4>q and </>i are localized in the ball 
Kr = {x G : |x| < R}, then by finite dependence on the domain of the signal in three di- 
mensions (Hyugen's principle), one obtains </>*(x, i) = for all x G Kr and t > t.R. Thus we have 
that 

(d t + Ud x )^*]=0, xGfi, t>t R . 

Thus it remains to consider flow variable </>**, whose aeroelastic potential on the boundary coincides 
with that of </>, and hence it displays regularity as in (|5.2p . This allows one to perform calculations as 

4 The general theory will provide at most H 2 / 3 (R? X [0,T]). 



in n X (0,T), 

on dn x (0,T), 
in R 3 + x (0,T), 

on % 9)) x(0,T). 
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in [HI Theorem 6.6.10]) in order to obtain the representation 



(d t + Ud x ft[</>] 



d(x,y,t) 



— ds d9[Mgd](x - (U + sin6)s,y - scos6,t - s) 
2tt Jo J 



1 f* 



+ 



Now using the same calculations as in [21 p. 333] we arrive at the following equation: 



u tt + A 2 u + ku t - [u, v{u) + F ] =Po- (d t + Ud x )u - q u (t) 



(5.3) 



for t large enough, with 




27r Jo Jo 



s,y — s cos 9, t — s). 



(5.4) 



Here, u is the extension of u by outside of il; Mg = sin9d x + cos9d y and t* is given by (|2.9[) . □ 

6 Appendix: Properties of delayed force q 

In this Appendix we consider properties of the delayed (aerodynamic type) force and prove Proposi- 
tion JO] and Lemma |2~T21 

6.1 Hidden Compactness of Retarded Potential: Proof of Proposition [2.41 

The proof of the bounds (f2~Kty - (f2~T2f can be found in [8] and [H]. Thus we need to check ([2~T3]) only. 
Without loss of generality we can assume u € C(— t*, +oo; C^°(il)). 
For any i/j e Hq(CI) we have 



recalling that x({7, 9, s) = (x—(U + sin 9)s, y — s cos 9). In all integrals above we extend the integration 




(6.1) 
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over Q to all of Mr and change spatial variables. 

,-27T 



[Mgu](T)tpdx dO 



[M$u](t - t*)ip(x(U,9, -t*))dx dO 



t" p2-K 



'0 JO 



<Qt(t),il>> 



(U + sm0)[M|u] x (x,r - s)ip(x(U,0, -s))dx d0 ds 

cos9[Mgu] y (x, t - s)ip(-x(U, 9, -s))dx d9 ds X. 

Now, we note that Mg can be moved under the integration in x, and we have |Mg/(x)| < \f x (x)\ + 
\f y (x)\ for all /. Hence, we integrate by parts with a single Me in both of the first integrals; in the 
third and fourth integrals, we also integrate by parts once as well. This leaves us with: 



V\ i- 



| <??(*), V> |<c{||ti(t)|| 1 + |Kt-t*)lli+ J , ll«(* + T-)lladr}| 
This implies the conclusion in (|2.13|) . The proof of Proposition 12.41 is complete. 
6.2 Proof of Lemma I27L21 



The relation in (|2.2ip easily follows from (|2.15|) and simple formula: 



dr (j){s)ds<t* <j>(s)ds, V <j> 6 Li(0,T). 

J T — t* J — 

Now we prove (|2.22p . Integrating by parts in £, and then applying (|2.17p with ip = u(t), we have: 
J <q u (T),u t (T)>dr\ <c{^ ||«(r)|| 2 ||u(r)|| a _„dT + jT \\u(r - t*)\\ 2 \\u(r)\\ 2 ^dT 

\\u(t - s)\\2\\u(r)\\ 2 -r,dsdr 



JO 



|g u (t)||-«r|K*)IU + ||g u (0)||_a||«(0)|| <r } 



<e ||u(r)||i + C e rsup||u(r)||l_ 



V 

-t* o,t 

\\u(t - s)\\jds dr + C e t* / \\u(t)\\ 2 2 dr 



o Jo 



n 



+ C{e\\q u (t)\\l a + - e \\u(t)\\l + e\\q»(0)\\-* + \\\u(0)\\l}. 
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After rescaling of e this implies 



<q u (r),u t (r)>dT < C(t*,e)(l + T) sup \\u{r)\\i_ v + e / \\u(p)\\ 2 2 dpdr 

rG[0,i] J-t* 
+ e{\\ q u (t)\f_ a + \\ q u (0)\\ 2 _ (T }. 

Therefore by (|2. 16|) this implies (|2.22j) with 77* = min{?7, 2 - a}. 
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